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abstract 

We present a new method for finding principal null directions (PNDs). Because our 
method assumes as input the intrinsic metric and extrinsic curvature of a spacelike hyper- 
surface, it should be particularly useful to numerical relativists. We illustrate our method 
by finding the PNDs of the Kastor-Traschen spacetimes, which contain arbitrarily many 
Q = M black holes in a de Sitter back-ground. 

Key Words: General Relativity, Asymptotic structure. Exact solutions. Numer- 
ical Relativity, Gravitational waves: theory. 



^Electronic address: shinkai@cfi.waseda.ac.jp 
^Electronic address: maeda@cfi.waseda.ac.jp 



1 Introduction 



According to General Relativity, light propagating from a spherical source to a distant 
observer through spacetime curved by the mass-energy of intervening bodies conveys to 
that observer a distorted (i.e., non-circular) image of the source. The amount of observed 
distortion depends on both the direction and the distance to the source, and of course 
vanishes as that distance approaches zero. Interestingly, the distortion growth rate also 
vanishes in this limit, independently of the direction from which the observer approaches 
the source. We may sum all this up by writing D = |C(e, 6, 0)e^, where D is a measure 
of the image distortion, e is a measure of the distance to the source, and the coefficient 
C(e, 6, (f)) depends smoothly on e (even at e = 0). Remarkably, for certain special directions 
{6i, (pi), it happens that C(0, 6, (p) = 0. In these directions, of which there can be at most 
four, the distortion vanishes unusually fast as the observer approaches the source. These 
approximately distortion-free directions are called principal null directions (PNDs) 

PNDs provide a great deal of gauge-invariant information about solutions to Einstein's 
equation. Indeed, in vacuum, they provide nearly as much information as the entire Rie- 
mann curvature. In the early 1960s, many researchers exploited this fact to obtain a 
large class of exact vacuum solutions with fewer than four distinct principal null direc- 
tions at each point P]. Some of these so-called algebraically special solutions were easy 
to interpret. Those with only one PND at each point, for example, clearly represented 
gravitational waves. But the significance of algebraically special solutions as a class re- 
mained somewhat obscure until Sachs proved his celebrated 'peeling theorem' 0]. This 
theorem asserts that the PNDs determined by the gravitational field of a bounded source 
gradually coalesce as distance from the source increases. The greater the distance, the 
more accurately we may approximate the field by ever more 'special' algebraically special 
solutions. 

Sachs's peeling theorem was only the first of many results linking the behavior of a 
spacetime's PNDs to the physical situation that the spacetime models. More recently, 
for example, Arianrhod et. al. have shown how to infer some interesting physical char- 
acteristics of the static, cylindrically symmetric Curzon spacetime from the (admittedly 
somewhat complex) behavior of its PNDs |Q. 

Despite their proven usefulness in constructing and interpreting exact solutions to 
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Einstein's equation, PNDs have played almost no role in work on numerical solutions. 
A language barrier seems at least partly responsible for this unfortunate state of affairs. 
Spinors and null tetrads have traditionally been the tools of choice in work on PNDs, and 
neither of these tools fits naturally into the '3+1' framework that is the mainstay of 
numerical relativity. 

In this paper, we present a new method for calculating PNDs tailored to the needs 
of numerical relativists. Our method assumes as input the induced metric and extrinsic 
curvature of a spacelike hypersurface, and produces as output the projections into that 
hypersurface of the PNDs. 

In the next section, we present the details of our method, and in section 3 we apply 
it to obtain the PNDs of the Kastor-Traschen solutions. The result is a series of pictures 
analogous to those in Arianrhod et. al. showing how the PNDs of various Kastor-Traschen 
solutions vary from point to point within the spacelike hypersurfaces of a natural slicing. 

2 A '3+1' Method for computing PNDs 
2.1 Our strategy 

Here, we present our method as it applies to vacuum spacetimes with cosmological con- 
stant A. We begin by fixing a triple (S, hab,Pab), where S is a smooth 3-manifold, hab is a 
Riemannian metric on S, and Pab is the extrinsic curvature of S. The constraints on hab 
and Pab are 

R-PabP^'+p' = 2A, (1) 
Daip'^'-ph'^') = 0. (2) 

Here p = Pabh"''', Da is the (unique) torsion-free derivative operator compatible with hab, 
and R = Rabh""^, where RabV^ = — 2D[aDm]f ™ for all smooth v"" on S. 

From {hab, Pab) we construct two further tensor fields Eab, Bab as follows: 

2 

Eab = Rab - Pa'^Pbm + PPab " "A/lafe, (3) 

Bab = er'DraP^^b, (4) 

where the tensor field Eabc = £[abc] satisfies Sabc^"''^'^ = 3! It follows from (0) and (0) that 
the fields Eab, Bab are both trace-free and symmetric. 
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The next step is to choose a unit vector field on S, and to decompose Eab, Bab into 
components along and perpendicular to z"". We set 

(5) 
(6) 

^eSab, (7) 
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(9) 
(10) 



where Sab = hab — ZaZ^- Finally, we set 



^0 = {-eab + Ja'bbc)m''m\ (11) 
^1 = ^(e.-J/6.)m^ (12) 

^2 = \{e-ib), (13) 
^3 = -^(e, + J/6,)m^ (14) 
^4 = {-Cab- Jabbc)m''m\ (15) 

where J a rotation by 90 degrees in the plane orthogonal to i", and = 

^(x" — iy"') for some pair of orthogonal unit vector fields that span that plane. 
Given \E'o5 ■ " " ? we need only solve the equation 

^4-2^ + 4^3^^ ^ 6^2^^ + 4^i;2 + $0 = 0. (16) 

to get our final answer: for each root Zi = tan ^e~^'^^ (i = 1, ■ ■ ■ , 4), the unit vector 

P(j)'' = cos Oiz"" + sin 9i cos ^jx" + sin 9i sin (17) 

determines a principal null direction. By this we mean the following: at points of S in 
the 4-dimensional maximal evolution Ai of (S, hab,Pab), t°' + P(i)°' is a principal null vector 
for each i, where is the unit normal to S in TW. 

This completes the description of our method. In the following section, we recall the 
steps required to solve 
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2.2 Procedure for finding PNDs 

We follow the d'lnverno-Russel- Clark method [|5[ to find the solutions Zi of (|T6|). Upon 
setting y = "if^z + \l/3, ([16D becomes 

+ QHy^ + AGy + K = 0, (18) 

where H, G and K stand for the following combinations of the complex numbers \l/o, ■ ■ ■ , ^1/4: 



/ 

J 

K 
H 

and G 



^4^0 - 4^i^3 + 3*2) 
\I^4 \&3 \I^2 

det *3 *2 *i 
*2 *i *o 



= *4* 



4^2 



*4*1 - 3*4*3*2 + 2*3. 



(19) 
(20) 

(21) 
(22) 
(23) 



The solutions of (|T8[) are simply expressed in terms of the solutions of 



X-^ -IX + 2J = 0, 



(24) 



which are 



Ai 
A2 
A3 



3-P 



(25) 
(26) 
(27) 



where P = {J + ^ .P - {I/?>f Y'^. From the A^ {i = 1,2,3) we now determine three 
further complex numbers, a, (3, 7, using the following equations: 



a 
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and a/37 



2*4Ai - AH, 
2*4A2 - 4:H, 
a^ + (3^ + AH, 
AG 



(28) 
(29) 
(30) 
(31) 



(p8|) - ( pOD determine a, [3, 7 up to sign, and (|31|) determines the signs. Then the solutions 
of (p!8|) are given as |(a + /5 + 7), |(a — /3 — 7), |(— a + — 7) and \{—0i — /? + 7). Finally, 



from a, f3, 7, we obtain the following four complex numbers: 

Zl = -{^3 + ^(a + /9 + 7)}/^4, (32) 

Z2 = -{vi>3 + ^(«-/3-7)}/^4, (33) 

= -{^3 + ^(-a + /5-7)}/^4, (34) 

Z4 = -{^3 + ^(-a-/3 + 7)}/^4. (35) 

These numbers are the solutions of (|T^). 

3 PNDs of the Kastor-Traschen solutions 

In this section, we illustrate our method by finding the PNDs of the Kastor-Traschen 
(KT) solutions |^. Because these solutions admit no timelike Killing vector field and 
contain black holes that undergo rapid relative motion, we would expect their PNDs to 
exhibit a variety of interesting and instructive behaviors. 

3.1 Kastor-Traschen solutions 

The KT solutions to Einstein's equation with cosmological constant contain arbitrary 
many Q = M black holes that participate in an overall de Sitter expansion or contraction. 
In the A ^ limit, the KT solutions reduce to the Majumdar-Papapetrou solutions, in 
which the balance between gravitational attraction and electrostatic repulsion among 
the black holes causes each to maintain its position relative to the others eternally. To 
write the KT metric, we first choose {xi,yi,Zi) G R'^, i = 1,2, ■ ■ ■ , N, and set rj = 
\/{x - Xi)2 + {y - + {z - Zif. Then 

ds^ = -^dt^ + aitWidx^ + dy"^ + dz^), (36) 

^ M- [a 
where = 1 + V — , a = e^* and H = ±\ -. 

V 3 

Naively, we interpret Mi as the mass of the iih black hole, although we have neither an 
asymptotically flat region nor event horizons available to convert this naive interpretation 
into a rigorous one. 
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3.2 PNDs of Kastor-Traschen solutions 

To apply the method of section 2 to the KT solutions, we must first choose spacelike 
hypersurfaces S and orthonormal triads x", y", z"- tangent to these S. The natural choice 
would seem to be the t = const, surfaces for E, and 

r = — f-V r = — (-Y z^ = — (- 

a^l \dx J ' aVl\dy)^ aVt\dz 

for the triad. We keep this choice in force throughout in what follows. 

We now select various values of the free parameters Zj) and Mj appearing in 

the KT metric (^), and plot the PNDs our method generates. For simplicity, we take all 
Zi = — that is, we confine all black holes to the x — y coordinate plane — and plot the 
projections of the PNDs into that plane. 

For the case of one black hole, the KT solution is just the Reissner-Nordstrom-de Sitter 
solution in cosmological coordinates, and its horizons are located at r-t = (1 — 2M\H\ ± 
^Jl — 4:M\H\) /2a\H\, where r_ and r+ denote the outer black hole horizon and de Sitter 
horizon, respectively. 

Figure 2 shows the PNDs of a KT solution containing two black holes of equal mass 
(Ml = M2 = 1.0) separated by a coordinate distance d = 5. In this case, r_|_ = 7.87M and 
r_ = 0.127M, so that the size of each black hole is the same as black points indicating 
its location. Four PNDs at each point are presented using arrows, but because these four 
PNDs coincide in pairs, with the coincident pairs pointing in opposite directions, we get 
a picture that looks like the electric field lines of equal charges. This indicates that the 
spacetime is Petrov type D at all points. 

As we decrease the separation between two black holes, we expect that at a certain 
critical value of that separation, their horizons will coalesce. When two equal-mass black 
holes inhabit a contracting {H < 0) KT spacetime (total mass M), the numerical results 
of p suggest a value of ~ 0.45M for this critical separation, which leads us to beheve 
that the two black holes in Figure 2 have distinct horizons. 

In Figure 3, we show the case of two black holes with different masses (Mi : M2 = 
1 : 10) separated by (a) d = 5, (h) d = 2 and (c) d = 0.5 in a contracting background 
{H < 0). Since the black holes move along the natural trajectories in that background, 
we may interpret these figures as snapshots of coalescing black holes. These cases have 
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none of the symmetry of Figure 2, and as a result, the field of PNDs is algebraically 
general (Petrov type I). However, we also see that the spacetime near each black hole 
is almost Petrov type D, as is a single black hole, and that at great distances from the 
holes the PNDs exhibit the same almost type-D behavior. This is just as we would expect. 
Moreover, we see that the PNDs on the axis connecting the two black holes exhibit Petrov 
type D behavior. Although we have no value for the critical separation in the case of two 
black holes with different masses. Figure 3(c) suggests that as the separation tends to 
zero, the spacetime's Petrov type approaches that of a single black hole. 

We can also draw the pattern of PNDs for KT solutions containing many black holes. 
As an example, we show in Figure 4 the PNDs of three black holes of equal mass located 
at the vertices of an equilateral triangle. 

4 Conclusions and discussion 

The results of the previous section illustrate the advantages of plotting PNDs over other 
methods of assessing how nearly algebraically special a spacetime is. Particularly when 
the spacetime is known only approximately, other methods are hard to apply, for they 
require us to determine whether certain scalar combinations of curvature components 
(e.g., — 6J^) are exactly zero or not. 

We hope the method presented in this paper will help numerical relativists to deal 
with at least two subtle issues. First, there is the problem of how to set initial data 
representing bounded sources free of incoming radiation. Because Sachs's peeling theorem 
assumes no incoming radiation, any failure of the PNDs computed from (E, hab,Pab) to 
exhibit the proper peeling behavior shows that (S, hab,Pab) is an unsatisfactory model of 
the instantaneous state of an isolated system. Our method thus allows us to recognize 
and reject initial data cluttered with spurious incoming radiation. This should be of some 
help. Second, there is the problem of detecting spurious gravitational waves in numerically 
generated spacetimes that result from physically inappropriate boundary conditions. We 
hope that such waves will manifest themselves in a pathological time dependence of the 
PNDs near numerical boundaries. 
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Figure Captions 



Figure 1: 

The solutions Zk (p^)-(p5|) for each k = 1, ...,4 at each point determine the PNDs. If we 
write Zk as = r^e"*''^'-', and set = tan then each PND P^^.^ points in the direction 
{0k, 0fc) indicated by the bold arrow in the figure, which we write as in (|1^. 

Figure 2: 

PNDs of the Kastor-Traschen solution for the case of two black holes of equal mass 
(Ml = M2 = 1) separated by a coordinate distance d = 5. Black holes are located at 
{x,y) = {3.0, 5.5) and (8.0, 5.5), and the arrows indicate PNDs at each point. 

Figure 3: 

PNDs of a KT solution containing two black holes with different masses (Mi : M2 = 1 : 10) 
separated by (a) = 5, (b) = 2 and (c) d = 0.5. In these cases, the symmetry of Figure 
2 is absent, we see that the field is algebraically general (Petrov type I). We also see, 
however, that the spacetime near each black hole and along the axis joining the two black 
holes is almost Petrov type D, and that in regions far from both holes, the spacetime is 
also almost Petrov type D, like a single black hole. 

Figure 4: 

PNDs of a KT solution with three black holes of equal mass at the vertices of an equilateral 
triangle. Black holes are located at {x,y) = {2.5, 3.0), (8.5, 3.0) and (5.5, 8.2). 
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